Acta Polytechnica Vol. YY No. X/ZZZZ 



The First Stars 

Daniel J. Whalen 1 

1 Me Williams Fellow, Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213 
Corresponding author: dwhalen@lanl.gov 

Abstract 

Pop III stars are the key to the character of primeval galaxies, the first heavy elements, the onset of cosmological 
reionization, and the seeds of supermassive black holes. Unfortunately, in spite of their increasing sophistication, numer- 
ical models of Pop III star formation cannot yet predict the masses of the first stars. Because they also lie at the edge 
of the observable universe, individual Pop III stars will remain beyond the reach of observatories for decades to come, 
and so their properties are unknown. However, it will soon be possible to constrain their masses by direct detection 
of their supernovae, and by reconciling their nucleosynthetic yields to the chemical abundances measured in ancient 
metal-poor stars in the Galactic halo, some of which may bear the ashes of the first stars. Here, I review the state of 
the art in numerical simulations of primordial stars and attempts to directly and indirectly constrain their properties 
observationally. 
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1 The Simulation Frontier 

Unlike with star formation in the Galaxy today, there 
is little disagreement over the initial conditions of star 
formation in the primordial universe. The original 
numerical simulations suggested that Pop III stars 
formed in small pregalactic structures known as cos- 
mological minihalos at z ~ 20 - 30, or ~ 200 Myr after 
the Big Bang (Bromm, Coppi & Larson, 1999,2001; 
Abel, Bryan & Norman, 2000, 2002; Nakamura & 
Umemura, 2001). These models predicted that the 
stars formed in isolation, one per halo, and that they 
were 100 - 500 M Q . Pop III stars profoundly trans- 
formed the halos that gave birth to them, expelling 
their baryons in supersonic ionized flows and later 
exploding as supernovae (e.g. Whalen, Abel & Nor- 
man, 2004; Kitayama k Yoshida 2005; Whalen & 
Norman 2008; Whalen et al. 2008a). Radiation 
fronts from these stars also engulfed nearby halos, 
either promoting or suppressing star formation in 
them, thereby regulating the rise of the first stellar 
populations (Shapiro et al. 2004; Iliev et al. 2005; 
Susa & Umemura 2006; Whalen et al., 2008b,2010; 
Hasegawa et al. 2009; Susa et al. 2009). 

The original estimates of Pop III stellar masses 
were not obtained by modeling the actual formation 
and evolution of the stars. They were derived by 
comparing infall rates at the center of the halo at 
early stages of collapse to Kelvin-Helmholz contrac- 
tion times to place upper limits on the final mass 
of the star. Later simulation campaigns in the same 
vein revealed a much broader range of final masses for 
Pop III stars, 30 - 300 M Q (O'Shea & Norman, 2007), 



and that they could form as binaries in a fraction of 
the halos (Turk et al., 2009). Heroic numerical efforts 
have only recently achieved the formation of a hydro- 
static protostar at the center of the halo (Yoshida et 
al., 2008) and the collapse of the central flow into an 
accretion disk (Stacy et al., 2010; Clark et al., 2011; 
Smith et al., 2011; Greif et al., 2011,2012). 



t= 5.56 yr 
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Figure 1: The formation and fragmentation of a Pop 
III protostellar disk in the Arepo code (Greif et al., 
2012). 

In particular, the disk calculations indicate that 
they are unstable to fragmentation, raising the pos- 
sibility that Pop III stars may have only been tens of 
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solar masses, not hundreds, and that they may have 
formed in small swarms of up to a dozen at the cen- 
ters of primeval halos. Computer models of ionizing 
UV breakout in the final stages of Pop III protostellar 
disks have also found that the I-front of the nascent 
star exits the disk in bipolar outflows that terminate 
accretion onto the star and mostly evaporate the disk 
by the time the star reaches ~ 40 M Q (Hosokawa et 
al., 2011). This result reinforces the sentiments of 
some in the community that while the Pop III IMF 
was top-heavy, primordial stars may only have been 
10 - 40 Ma. 



1.1 High-Mass or Low-Mass Pop III 
Stars? 

In spite of their increasing sophistication, these sim- 
ulations should be taken to be very preliminary for 
several reasons. First, Pop III accretion disks form 
in smoothed-particle hydrodynamics (SPH) models 
but not in adaptive mesh refinement (AMR) simu- 
lations, although the AMR models have not evolved 
the collapse of the halo to the times achieved by SPH 
calculations. This raises the question of whether one 
technique better captures the transport of angular 
momentum out of the center of the cloud than the 
other, and whether accretion is ultimately spherical 
or through a disk. Second, the stability of the disk 
itself remains an open question because although the 
simulations can now fully resolve the disk they do not 
yet incorporate all of its relevant physics. In partic- 
ular, they lack high-order radiation transport, which 
regulates the thermal state of the disk and its ten- 
dency to fragment. Furthermore, the role of primor- 
dial magnetic fields in the formation and evolution 
of the disk is not well understood (Turk et al., 2012; 
Widrow et al., 2012). 

The use of sink particles to represent disk frag- 
ments in the original SPH simulations of Pop III 
protostellar disk formation called into question the 
longevity of the fragments. Once they are created in 
the simulations they are never destroyed, unlike real 
fragments that could be torn apart by gravitational 
torques and viscous forces (Norman, 2010). More 
recent moving mesh simulations performed with the 
Arepo code that do not rely on sink particles find that 
the fragments persist but only evolve the disk for 10 
- 20 yr (Greif et al., 2012). Perhaps most impor- 
tantly, no simulation has followed the disks for more 
than a few centuries, far short of the time required 
to assemble a massive star. Thus, it remains unclear 
if the fragments in the disk remain distinct objects 
or merge with the largest one at the center, build- 
ing it up into a very massive star over time through 
protracted, clumpy accretion. 



1.2 Accretion Cutoff and the Final 
Masses of Pop III Stars 

This latter point directly impacts estimates of final 
masses for Pop III stars inferred from numerical sim- 
ulations that attempt to model how ionizing UV from 
the star reverses infall and evaporates the accretion 
disk. At the heart of such models is a simple recipe 
for the evolution of the protostar that provides a pre- 
scription for its radius and luminosity as a function of 
time and acquired mass. The Hosokawa et al. (2011) 
2D calculations take the growth of the protostar to 
be relatively steady, in which case it contracts and 
settles onto the main sequence at ~ 30 M . At this 
point the star becomes extremely luminous in ioniz- 
ing UV radiation that halts accretion onto the star in 
a few hundred kyr at a final mass of ~ 40 M . If ac- 
cretion instead turns out to be clumpy, the protostar 
could remain puffy and cool and reach much larger 
masses before burning off the disk. 

A finer point is that all current accretion cutoff 
simulations evolve both radiation and hydrodynam- 
ics on the Courant time, a practice which is known 
to lead to serious inaccuracies in I-front propaga- 
tion in density gradients (Whalen & Norman, 2006). 
Such coarse time steps may result in premature I- 
front breakout and accretion cutoff, and hence un- 
derestimates of the final mass of the star. Three di- 
mensional simulations with more accurate radiation- 
matter coupling schemes, both steady and clumpy 
accretion scenarios, more realistic prescriptions for 
protostellar evolution based on nucleosynthesis codes 
such as KEPLER (Weaver et al., 1978; Woosley et al., 
2002) and a variety of halo environments may better 
constrain the Pop III IMF. However, in judging the 
power of such simulations to model the masses of the 
first stars, it should be remembered that no simula- 
tions realistically bridge the gap in time between the 
formation and fragmentation of a protostellar disk 
and its photoevaporation up to a Myr later. We note 
in passing that fragments can also stop accreting if 
they are ejected from the disk by 3-body gravita- 
tional effects (Greif et al., 2011; Johnson & Khochfar, 
2011). These fragments could become very low-mass 
Pop III stars (~ 1 M Q ); if so, some of them may live 
today. 

2 Constraining the Pop III 
IMF with Stellar 
Archaeology 

Unfortunately, because they lie at the edge of the 
observable universe, individual Pop III stars will 
remain beyond the reach of direct detection for 
decades to come, even with their enormous lumi- 
nosities (Schaerer 2002) and the advent of the next 



Acta Polytechnica Vol. YY No. X/ZZZZ 




Figure 2: Comparing Pop III SN yields to the chemical abundances of three of the most metal-poor stars (left 
panel) and the extremely metal-poor (EMP) stars in the Cayrel et al. (2004) and Lai et al. (2008) surveys 
(right panel). In the left panel, the abundances in HE0557-4840 agree well with the yields from SN model zl5G 
in the Joggerst et al. (2010) study. In the right panel we show that higher explosion energy rotating Z = 
stars reproduce EMP abundances well. The existence of 15 M© Pop III stars is required to produce this good 
agreement with observations. 



generation of near infrared (NIR) observatories such 
as the James Webb Space Telescope (JWST) and 
the Thirty-Meter Telescope (TMT). However, there 
have been attempts to indirectly constrain the masses 
of Pop III stars by comparing the cumulative ele- 
mental yield of their supernovae to the fossil chemi- 
cal abundances found in ancient metal-poor stars in 
the Galactic halo, some of which may be second- 
generation stars. Stellar evolution models indicate 
that 15 - 40 Mq primordial stars died in core collapse 
(CC) supernovae (SNe) and that 40 - 140 M© stars 
collapsed to black holes, perhaps with violent pul- 
sational mass loss prior to death (Heger & Woosley, 
2002). Pop III stars between 140 and 260 M can 
die in pair-instability (PI) SNe, with energies up to 
100 times those of Type la SNe that completely un- 
bind the star and leave no compact remnant (Chat- 
zopoulos & Wheeler 2012 have recently discovered 
that rotating Pop III stars down to 65 M Q can also 
die as PI SNe) . These explosions were the first great 
nucleosynthetic engines of the universe, expelling up 
to half the mass of the progenitor in heavy elements 
into the early IGM. Primordial stars above 260 M© 
collapsed directly to black holes, with no mass loss. 

Joggerst et al. (2010) recently calculated the 
chemical imprint of low-mass Pop III SNe on later 
generations of stars by modeling mixing and fallback 
onto the central black hole in 15 - 40 M© Pop III 
core collapse explosions with the CASTRO AMR 
code. As shown in Figure [2j a simple power-law 
IMF average of the elemental yields of these explo- 



sions is in good agreement with the fossil abundances 
in a sample of 130 extremely metal poor stars with 
Z < 10" 4 Z Q (Cayrel et al., 2004; Lai et al., 2008). 
Although these results suggest that low-mass Pop III 
stars shouldered the bulk of the chemical enrichment 
of the early IGM, 40 - 60 M© hypernova explosions, 
whose energies are intermediate to those of CC and 
PI SNe, may also have contributed metals at high 
redshifts (Iwamoto et al. 2005). 

To date, the telltale odd-even nucleosynthetic sig- 
nature of PI SNe has not been found in the fossil 
abundance record, leading some to assert that Pop 
III stars could not have been very massive. How- 
ever, the odd-even effect may have been masked by 
observational bias in previous surveys (Karlsson et 
al., 2008). Reconciling Pop III SN yields to the ele- 
mental patterns in metal-poor stars is still in its in- 
fancy for several reasons. First, only small numbers 
of extremely metal-poor stars have been discovered to 
date, and larger sample sizes would better constrain 
early SN yields. Second, measurements of some el- 
ements in low-metallicity stars are challenging and 
in the past have been subject to systematic error. 
Finally, there are many intervening hydrodynamical 
processes between the expulsion of the first metals 
and their uptake into new stars that are not yet un- 
derstood. 
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Figure 3: Pop III PI SN light curves for the JWST NIRCam. Clockwise from the upper left panel, the redshifts 
are z = 10, 15, 20 and 30. The optimum filter for each redshift is noted on the y-axis labels and the times on 
the x-axes are for the observer frame. F200, F277 and F356 are at 2.0, 2.77 and 3.56 /zm, respectively. 



3 Finding the First Cosmic 
Explosions 

Detection of Pop III SNe would unambiguously probe 
the masses of primordial stars for the first time. Since 
these explosions are 100,000 times brighter than ei- 
ther their progenitors or the primitive galaxies that 
host them, they could be found by JWST or the 
Wide-Field Infrared Survey Telescope (WFIRST). 
However, unlike the Type la SNe used to constrain 
cosmic acceleration, light from primeval supernovae 
must traverse the vast cosmic web of neutral hydro- 
gen that filled the universe prior to the epoch of reion- 
ization. Lyman absorption by this hydrogen removes 
or scatters most of the light from ancient supernovae 
out of our line of sight, obscuring them. 

Whalen et al. (2012a,b,c) have calculated JWST 
NIR light curves for Pop III pair-instability SNe with 



the Los Alamos RAGE and SPECTRUM codes (Git- 
tings et al., 2008; Frey et al., 2012), which are shown 
for z — 10, 15, 20 and 30 in Figure [3] These simula- 
tions include radiation hydrodynamical calculations 
of the SN light curve and spectra in the local frame, 
cosmological redshifting, and Lyman absorption by 
intergalactic hydrogen. JWST detection limits at 2 
- 4 /xm are AB magnitude 31 - 32, so it is clear that 
JWST will be able to detect the first cosmic explo- 
sions in the universe if they are PI SNe (and even 
perform spectrometry on them). Even given JWST's 
very narrow fields of view at high redshifts, recent cal- 
culations indicate that at least a few PI SNe should 
be present in any given JWST survey (Hummel et al. 
2012). Also, because WFIRST detection limits will 
be AB magnitude 26.5 at 2.2 ^m, it is clear from Fig- 
ure H that Pop III PI SNe will be visible to WFIRST 
out to z ~ 15 - 20. Since it is an all-sky survey, and 
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because this redshift range may favor the formation 
of very massive Pop III stars because of the rise of 
Lyman- Werner UV backgrounds (O'Shea & Norman, 
2008), WFIRST will detect much larger numbers of 
Pop III SNe. 

Could Pop III SNe be detected at later stages by 
other means? Whalen et al., (2008b) found that most 
of the energy of Pop III CC SNe is eventually radiated 
away as H and He lines as the remnant sweeps up and 
shocks surrounding gas. At later epochs this energy 
would instead be lost to fine structure cooling by met- 
als. In both cases the emission is too dim, redshifted 
and drawn out over time to be detected by any up- 
coming instruments. However, PI SNe deposit up to 
half of their energy into the cosmic microwave back- 
ground (CMB) by inverse Compton scattering at z ~ 
20 (Kitayama & Yoshida, 2005; Whalen et al., 2008b) 
and could impose excess power on the CMB at small 
scales (Oh, Cooray & Kamionkowski, 2003). The res- 
olution of current ground-based CMB telescopes such 
as the Atacama Cosmology Telescope and South Pole 
Telescope approaches that required to directly image 
Sunyaev-Zeldovich (SZ) fluctuations from individual 
Pop III PI SN remnants, so future observatories may 
detect them. Unlike PI SNe, CC SNe deposit little 
of their energy into the CMB at z ~ 20, and even 
less at lower redshifts because the density of CMB 
photons falls with cosmological expansion. 

The extreme NIR luminosities of primordial PI 
SNe could contribute to a NIR background excess, 
as has been suggested for Pop III stars themselves, 
i.e. Kashlinsky (2005). New calculations reveal that 
enough synchrotron emission from CC SN remnants 
is redshifted into the 21 cm band above z ~ 10 to 
be directly detected by the Square Kilometer Ar- 
ray {SKA) (Meiksen & Whalen 2012). Somewhat 
more energetic hypernovae could be detected by ex- 
isting facilities such as the Extended Very-Large Ar- 
ray (eVLA) and eMERLIN. PI SN remnants gener- 
ally expand into ambient media that are too diffuse to 
generate a detectable synchrotron signal. Pop III SN 
event rates make it unlikely that they will be found 
in absorption at 21 cm at z > 10. 

The detection of the first cosmic explosions will 
be one of the most spectacular discoveries in extra- 
galactic astronomy in the coming decade, opening our 
first observational window on the era of first light 
and the end of the cosmic Dark Ages at z ~ 30. 
They will unveil the nature of primordial stars and 
constrain scenarios for early cosmological reioniza- 
tion, the process whereby the universe was gradually 
transformed from a cold, dark, featureless void into 
the vast, hot, ionized expanse of galaxies we observe 
today. At somewhat lower redshifts (z ~ 10 - 15), 
detections of Pop III supernovae will probe the era 
of primitive galaxy formation, marking the positions 
of nascent galaxies on the sky that might otherwise 



have eluded detection by JWST. Finally, finding the 
first supernovae could also reveal the masses of the 
seeds of the supermassive black holes lurking at the 
centers of massive galaxies today. 
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DISCUSSION 

Wolfgang Kundt: You spoke of a mass range (> 40 
M Q ) for primordial stars for which they would col- 
lapse to black holes. Where are the nearest of them 
today? 

Daniel Whalen: This is an excellent question. We 
recently published a letter (Whalen & Fryer 2012) 
in which we found that most 20 - 40 M Q Pop III 
black holes would be ejected from the cosmological 
halos that gave birth to them at velocities of 500 - 
1000 km s _1 by natal kicks due to asymmetries in the 
core-collapse engine. Such velocities are far above the 
escape velocity of any halo they would encounter for 
over a Hubble time, so there is a good chance that 
these black holes would be exiled to the voids be- 
tween galaxies today. Black holes above 40 M Q are 
unlikely to be born with kicks and remain in the halo, 
accreting and growing over cosmic time. These black 
holes are much more likely to reside in the galaxies 
into which their host halos were taken, a few of which 
could become the supermassive black holes found in 
the SDSS quasars today. 

Maurice Van Putten: What fraction of the gas in 
a cosmological halo ends up in primordial stars? 

Daniel Whalen: It is currently thought that the 
minimum halo mass for forming a Pop III star is ~ 
10 5 M Q and that 1-10 stars are formed with masses 
of 30 - 300 M.q. Thus, a conservative estimate is that 
0.1 - 1% of the baryons in the halo are converted into 
stars, and that the rest are evicted from the halo by 
strong ionized flows over the life of the stars. 



